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We have calculated and statistically analyzed the magnetic-field spectrum 
(the "B-spectrum" ) at fixed electron Fermi energy for two quantum dot systems 
with classically chaotic shape. This is a new problem which arises naturally in 
transport measurements where the incoming electron has a fixed energy while 
one tunes the magnetic field to obtain resonance conductance patterns. The 
"B-spectrum", defined as the collection of values {B{\ at which conductance 
g(Bi) takes extremal values, is determined by a quadratic eigenvalue equation, in 
distinct difference to the usual linear eigenvalue problem satisfied by the energy 
levels. We found that the lower part of the "B-spectrum" satisfies the distribution 
belonging to Gaussian Unitary Ensemble, while the higher part obeys a Poisson- 
like behavior. We also found that the "B-spectrum" fluctuations of the chaotic 
system are consistent with the results we obtained from random matrices. 
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Due to recent advances in controlled crystal growth and lithographic tech- 
niques, it is now possible to fabricate various "artificial atoms" or quantum dots 
whose size is so small such that transport is in the ballistic regime 0. Among 
the many interesting phenomena associated with quantum ballistic transport, it 
was proposed that these quantum dots could be used to examine the theoret- 
ical notion of "quantum chaos" |3], [|, |5| . More recently, Taylor et. al. fabricated 
an electronic Sinai billiard and investigated fluctuations of the resistance as 
an external magnetic field is tuned. An important discovery of this work is the 
apparent fractal resistance fluctuations |7J which maybe discussed on the basis of 
a semiclassical theory ||. 

The problem of quantum chaos is a very interesting theoretical physics issue, 
and our theoretical understanding of it has been advanced by Random Matrix 
Theory (RMT) which classifies the statistics of the eigen-spectrum of a quantum 
system according to the Wigner-Dyson ensembles |J. It is well established that, 
for a classically non-chaotic system such as a particle confined to move inside 
a rectangular box, the normalized nearest neighbor energy spacings {s} obey a 
Poisson distribution [|IIJ] P(s) = e~ s . On the other hand, for a classically chaotic 
system such as a spinless particle confined to move inside a stadium shaped 



box[ll, 12], the spacings follow the Wigner distribution which belongs to the 



Gaussian Orthogonal Ensemble (GOE) in the language of RMTfllBf. Furthermore, 
when time reversal invariance is broken, say by applying a magnetic field, the 
system is described [[L4| by Gaussian Unitary Ensemble (GUE), where P-z{s) ~ 



32s" - 

-9. C 



In order to study certain aspects of quantum chaos using quantum transport 
techniques^, [15|, 0], one must deal with open systems where a scattering problem 
of charge carriers by some peculiar boundary must be solved. Experimentally it 
is quite difficult to study conduction as a function of the incoming electron en- 
ergy in a quantitatively accurate fashion |16]], although measurements on tunnel 
junctions made of Al nanoparticles have recently be made|T7| and its relation to 
quantum chaos discussed Using 2D electron gas fabricated with compound 
semiconductors, experiments usually measure conductance as a function of exter- 
nal magnetic fieldsp], 15| |6], g(B), at a fixed electron Fermi energy E D . When 



a quantum dot is weakly coupled to the external leads, the magneto-conductance 
g{B) may show resonance like behavior as the magnetic field B is varied||, if the 
measurement is indeed in a regime which probes the internal electronic states fll7|. 
This behavior may be understood as follows. As B is varied, the energy levels {e{\ 
of the scattering states labeled by indices i = 1, 2, ... in the quantum dot change 
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with it: €i = €i(B). These levels are well separated since the the dot region is 
weakly coupled to the leads. For an incoming electron with a fixed Fermi energy 
E Q , each time when the internal state energy e,(5) is tuned to equal to E Q as B 
is varied, a resonance peak occurs in g(B) due to a junction resonance. This was 
indeed observed in numerical simulations |J by solving the quantum scattering 
problem. Clearly this behavior should be observable if the Coulomb blockade 
effects are small, which happens when the system has a large capacitance thus 
small charging energy. We shall assume this to be the case. 

Hence in this junction resonance regime, it is interesting to define a "B- 
spectrum" as the collection of values of the special magnetic fields B = {Bi} 
at which g(Bi) is peaked f|. It is important to ask the following questions: what 
are the statistical properties of this "B-spectrum" ? How do its statistical prop- 
erties change with the system shape ? These are also useful questions to answer 
because it is increasingly possible to directly probe the internal electronic states 
of an isolated quantum dot, as demonstrated by the experiment reported in Ref. 

Motivated by these questions, in this paper we report our studies on closed 
quantum dot systems where these questions can be answered clearly. Simons, 
Szafer and AltshulerQ have investigated the correlations of slopes of the energy 
levels as a function of an external parameter such as the magnetic field B. We 
however emphasize that in their studies, the focus is on how a particular energy 
level Ei varies with the external parameter B. Our investigation, on the other 
hand, is completely different as it does not focus on any energy level: the energy 
E is a fixed parameter in the Schrodinger equation by the equilibrium Fermi 
energy of the system. Rather, we investigate the statistical properties of the set 
{Bi} which makes E a an eigenvalue of the Hamiltonian. Indeed, we note that 
fixing B and study energy levels {Ei}, or fixing E Q and study the "B-spectrum" 
{Bi}, are two different problems. The former is about the energy eigenvalue 
spectra and its relation to an external parameter B, while the latter is relevant for 
transport situations where the incoming electron has energy E which is fixed by 
the electron reservoir and cannot change, while one tunes B to special values {Bi} 
where conductance g(Bi) is peaked. To our knowledge the statistical ensemble 
which is satisfied by the "B-spectrum" has not previously been determined. As 
we shall see below, this "B-spectrum" is determined by a quadratic eigenvalue 
problem, in distinct contrast to the usual energy spectrum at a fixed B which is 
a linear eigenvalue problem. When a quantum dot is weakly coupled with the 
external leads, our answers are valid since in the weak coupling regime the levels 
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are well separated and statistical properties should not change |20], j2jj. 

To make the problem at hand clearer, the inset of Fig. (3) shows the two- 
dimensional quantum dot systems we have studied: a Sinai-like billiard and a 
stadium-shaped quantum dot. These systems are classically chaotic systems. 
Experimentally transport measurements have been reported through these quan- 
tum dots which were fabricated using the split-gate technology [0, [| |2|. Fig. 
(1) shows the energy levels of a Sinai-like billiard as a function of magnetic field 
B. For a given B = B Q , the levels satisfy GUE as mentioned above. These are 
the intersections of the line B = B Q (the vertical solid line) with the spectrum. 
However, for a fixed energy E = E , we are interested in the intersections of the 
horizontal solid line with the spectrum which define the "B-spectrum" . From the 
curves of Fig. (1) it is clearly not obvious what statistics the "B-spectrum" will 
satisfy. 

In the presence of a magnetic field, the single-particle Schrodinger equation 
can be written as 

{-V 2 -e + b 2 A 2 + 2biA - V}^ = , (1) 

where A = BAq = Bxj is the vector potential, B is the magnetic field, E = 
h 2 e /2m is the energy, b = eB/hc, and i 2 = —1. Discretizing the spatial deriva- 
tive, Eq. (HD can be cast into a matrix form, 

(M l + ibM 2 + b 2 M 3 )^ = e *lj , (2) 

where Mi, M 2 , M 3 are real matrices and M 3 is also diagonal. From Eq.(||D it 
is clear that for a fixed magnetic field the solution of all the allowed energies 
forms the usual linear eigenvalue problem which has been studied intensively. 
However for a fixed energy eo, the solution of all the allowed magnetic fields B, 
the "B-spectrum" , forms a quadratic eigenvalue problem. 

There are two methods to find the "B-spectrum". The first method is by 
brute force: one directly calculates the energy eigenvalues E for a given magnetic 
field B and traces out the curves of E versus B. These curves may cross over each 
other. One then calculates, from this E versus B curve, the set of magnetic fields 
B for a fixed energy E a . Fig. (1) shows the E versus B curves for the first 40 
eigenstates obtained using the Lanczos eigenvalue technique for the Sinai-like 



billiard system. However, this method gets time consuming very quickly if higher 
and higher states are needed. 

The second method is to transform the quadratic eigenvalue problem Eq. (^) 
into a usual linear eigenvalue problem [p3fl . Let ^(t) = e M ip and define N\ = 
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M 3 _1 (Mi - e ), N 2 = Mz 1 M 2 , so that Eq. (§) becomes + JV 2 tf' + = 

where = ib^ is the first derivative with respect to the parameter t, and / 
is the unit matrix. With vp' = $ ($' = ib<f>), we have $ = M>' = and 
iV 2 $ + = W = ib®, or 

T \ / iTr \ / lTr \ 

(3) 

This is a linear eigenvalue problem for b. The matrix in Eq. (|3|) is however not 
Hermitian, hence its eigenvalues may or may not be real. To obtain the physical 
solution, we look for all imaginary eigenvalues so that b is real. We verified that 
the results coincide with those obtained with the conventional Lanczos technique 
described above. In this way we have calculated the "B-spectrum" from Eq. ([3D 
for the two chaotic systems (inset of Fig. (3)), for different energies cq. The 
confining potential is assumed to be hard wall. 

For the Sinai-like billiard, we considered a 4.7/zmx4.3/im rectangular quantum 
dot with three hard disks inside the dot. The radius of the disks were fixed at 
0.9/im, 0.8/im, and 0.5//m with their centers randomly chosen but without the 
disks overlapping each other. Randomly changing positions of the disks allows us 
to generate different configurations for ensemble averaging to obtain reasonable 
statistics. The number of physical solutions N obtained from Eq.@ depends 
on the fixed energy eo, with N increasing with the value of the energy e®. For 
instance, when eo is fixed at 30meV, we have N ~ 1433 physical solutions for 
the "B-spectrum" out of a total of about 3000 eigenvalues. We found that the 
statistics of the "B-spectrum" behaves differently for the lower and higher part of 
the spectrum respectively. The lowest f24}| ~ 300 "B-levels" give GUE statistics 
and the highest ~ 800 "B-levels" give a Poisson-like behavior. Fig. (2) plots the 
distribution function obtained from our numerical data of the nearest neighbor 
"B-level" spacings for the Sinai-like billiards. With the ensemble average of 20 
different configurations, the distribution determined from the lower part of the 
"B-spectrum" agrees well with GUE statistics (solid line) ^(s) discussed above. 
On the other hand, the higher part of the spectrum, shown in the inset of Fig. 
(2), has a Poisson-like behavior. 

Another often used measure in studying level statistics is the spectral rigidity 
A 3 , defined as|| ^] the mean square deviation of the best local fit straight 
line to the staircase cumulative spectral density over a normalized energy scale. 
This quantity measures longer range correlations of the level spectrum and often 
provides a more critical test of the level statistics. To compute A 3 we followed 



a scheme presented in Ref. |[26|| , and the numerical data is compared with the 
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analytical formula from random matrix theory p7j. Fig (3) shows the A3 analysis 
of the lower part of the "B-spectrum" . It is clear that the data is in very good 



agreement with the GUE statistics 1 28 



To test the statistical properties of the "B-spectrum" further, we studied 
another chaotic system, namely a stadium-shaped quantum dot (inset of Fig. 
(3)). The distribution function and A3 for its "B-spectrum" are included in 
Figs. (2) and (3). Here we did not use an ensemble average and the data is for 
one system only. The general trend is the same as for the Sinai-like billiards, 
and still clearly shows the two distinct behaviors for different parts of the "B- 
spectrum" , namely a GUE behavior of the lower part and a Poisson-like behavior 
for the higher part. Finally, we have verified that the same statistical behavior 
is observed for both the Sinai-like billiard and the stadium-shape billiard with 
many different values of the fixed energy eo, and conclude that the lower part of 
the "B-spectrum" of both systems satisfies GUE statistics. 

It is not difficult to understand that the higher part of the "B-spectrum" 
should behave differently. The higher part corresponds to larger values of the 
magnetic fields, which are known to destroy chaos [p9| , [J. For the particular 
sizes of the Sinai-like billiard and the stadium-shaped dot, the magnetic-field B c 
that roughly separates the low lying and high lying part of the "B-spectrum" 
is about B c ~ 3 to 4T. The classical cyclotron radius of the electron at and 
above this field strength is quite small compared with the system size. Hence the 
electron "skips" along the wall of the confining potential or makes circular motion 
inside the quantum dot, thereby reducing the effect of chaotic scattering by the 
geometry [Q. Nevertheless, the field range up to B c is quite wide and it should be 
possible to investigate the "B-spectrum" experimentally using resonant magneto- 
conductance measurements for systems having weak coupling to the leads. 

While the quadratic eigenvalue problems reported above were for billiard sys- 
tems where the continuum Schrodinger equation is discretized to obtain the ma- 
trix equation (0), we also studied a similar quadratic eigenvalue problem using 
random matrices to replace the matrices Mi, M 2 and M3 in Eq. (||]). Two re- 
quirements must be satisfied: first, the matrices Mi, M 2 , and M 3 must be real; 
second, Mi must be symmetric and M 2 antisymmetric so that the Hamiltonian 
be Hermitian. The random matrices were set up in the standard fashion: for a 
N x N matrix, Mi has N(N + l)/2 independent matrix elements and M 2 has 
N(N — l)/2, where the matrix elements are Gaussian random numbers. We then 
diagonalized Eq.(|3]) using the same numerical methods discussed above and found 
all the physical solutions. We chose N = 1500 so that the matrix which must 
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be diagonalized is 3000 x 3000. Typically we obtained about 700 physical "B- 
levels" for various energies e . The results are included in Figs. (2) and (3). It 
is clear that the "B-level" statistics from the random matrices is completely con- 
sistent with the GUE statistics. We may conclude that "B-spectra" coming from 
the quantum billiard systems and from random matrices are still consistent with 
each other and a universality can still be established for this quadratic eigenvalue 
problem using the random matrices. 

In summary, we have numerically investigated the statistical properties of the 
magnetic field spectra (the "B-spectrum"), which is determined by a quadratic 
eigenvalue problem. This spectrum is defined by the allowed magnetic fields for 
an electron moving in a quantum dot with its fixed fermi energy. This prob- 
lem arises for systems with weak coupling to the leads in which the scattering 
states are well separated in energy. For two different chaotic billiards, e.g. two- 
dimensional quantum dots in the shape of a Sinai-like billiard and a stadium 
billiard, the "B-spectra" have distinctly different statistical behavior at the lower 
and higher parts of the spectra. In particular, the lower part is well described by 
the GUE statistics while the higher part is Poisson like. We found that the same 
quadratic eigenvalue problem can be studied using random matrices as well, and 
the eigenvalues from the random matrices have precisely the same behavior as 
those of the billiards. Thus the notion of universality classes using the random 
matrix theory can be carried over to this new problem. While our numerical 
data provided clear evidence of the statistical properties of the present quadratic 
eigenvalue problem, it is not at all obvious a priori that such statistical proper- 
ties are controlled by the GUE universality. Further work is needed to provide 
an analytical understanding. 

Experimentally the statistical properties of the "B-spectra" could be examined 
for quantum dot systems which couple weakly to the external leads and thus the 
transmission is controlled by junction resonances. The weak coupling could be 
provided by adding constrictions at the connections of the leads with the quantum 
dot for which only the lowest few subbands of the leads can propagate. In typical 
experimental situation on submicron structures the single particle level spacing 
is around 0.05meV, thus they can be measured if the temperature is kept to less 
than 500mK. We thus conclude that for a magneto-transport measurement in 
the junction resonance regime, the special magnetic field strengths {Bi} at which 
g(Bi) takes extremal values, should satisfy GUE statistics. 
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Figure Captions 

Figure 1. Typical curves of the energy levels as a function of the magnetic field, 
E = E(B), for the Sinai-like billiard. The vertical solid line at B = B Q 
intersects the spectrum and the intersections give the usual energy levels 
E = E(B Q ). The "B-spectrum" is the collection of intersections of the 
horizontal solid line at a fixed energy E = E Q . 

Figure 2. The distribution function of the nearest neighbor "B-level" spacings ob- 
tained from the lower part of the "B-spectrum". The solid line is the an- 
alytical formula for GUE ^(s)- Solid squares: ensemble averaged data 
for 20 independent Sinai-like billiards; open circles: for a stadium-shaped 
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quantum dot; triangles: for random matrices averaged over 13 configura- 
tions. Insert: the distribution function obtained from the higher part of 
the "B-spectrum" of the same systems (symbols) and the solid line is the 
Poisson distribution P(s). The energy is fixed at e = 30meV. 

Figure 3. The spectral rigidity A^(L), where L is the number of the "B- levels" in- 
volved in computing A3. The solid line is the analytical expression for GUE 
statistics, taken from Ref.[18]. Solid squares: ensemble averaged data for 20 
independent Sinai-like billiards; open circles: for a stadium-shaped quantum 
dot; triangles: from the random matrices averaged over 13 configurations. 
Insert: sketch of the chaotic structures studied here, a Sinai-like billiard 
with three hard disks confined inside a rectangle; and a stadium-shaped 
quantum dot. 
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